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Motivated by thermodynamic considerations, we analyse the variation of the quantum mutual 
information on a unitary orbit of a bipartite system's state, with and without global constraints 
such as energy conservation. We solve the full optimisation problem for the smallest system of 
two qubits, and explore thoroughly the effect of unitary operations on the space of reduced-state 
spectra. We then provide applications of these ideas to physical processes within closed quantum 
systems, such as a generalized collision model approach to thermal equilibrium and a global Maxwell 
demon playing tricks on local observers. For higher dimensions, the maximization of correlations 
is relatively straightforward for equal-sized subsystems, however their minimisation displays non- 
trivial structures. We characterise a set of separable states in which the minimally correlated state 
resides: a collection of classically correlated states admitting a particular "Young tableau" form. 
Furthermore, a partial order exists on this set with respect to individual marginal entropies, and 
the presence of a "see-saw effect" for these entropies forces a finer analysis to determine the optimal 
tableau. 
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I. INTRODUCTION 

The idealisation that a quantum system evolves in a 
closed environment dates back to the foundations of the 
theory itself [1]. While this closure preserves the pu- 
rity of any quantum state, the correlations between its 
constituent parts can greatly vary. A natural but little- 
addressed question is how much can these correlations 
vary when the system undergoes isolated, unitary dy- 
namics? The issue of varying correlations arises natu- 
rally in different thermodynamic processes within closed 
quantum systems, and can often result in seemingly para- 
doxical effects. 

We address this issue for a bipartite system and use 
the quantum mutual information (QMI) as a measure of 
total correlations. The QMI is a quantity that appears 
within a range of scenarios in quantum thermodynam- 
ics, such as work extraction, heat transfer and the ap- 
proach to equilibrium [2] . This paper builds on the work 
in Ref. [2] and illustrates that even the simplest bipar- 
tite systems display a rich and non-trivial structure that 
deserves careful study. 

We start in Section |ll] with a declaration of the core 
mathematical objects and notions that will be used 
throughout the paper. Section |ITT] formally states the 
key underlying task that we tackle, namely the variation 
of the quantum mutual information on a unitary orbit. 
We present the general solution for the maximum QMI 
in section |III A| and detail what is required for its min- 
imisation in section |IIIB| Section |III C| provides a full 
and explicit solution of this optimisation problem for the 
instance of two qubits, including the case when a global 
constraint (e.g. energy conservation) is present. An im- 
portant tool is introduced: the set TZ of eigenvalues for 
marginals compatible with all states on a unitary orbit of 
a bipartite quantum system. In section [Tv| we study this 
set TZ in more detail and observe the effect that unitaries 



(acting on the global state) have on the marginal spectra. 

With these results established, in sect ion [V| we describe 
some applications to physical processes within a closed 
quantum system. Section |V A| discusses the heat flow 
model of Ref. [2] in more detail and covers some im- 
portant assumptions. In section V C we establish a more 



general collision model for equilibration in which the sys- 
tems retain correlations. Finally, section [V D| illustrates 
how a global Maxwell-type demon can confuse a local 
observer by causing an apparent violation of his Second 
Law. 

Section [Vl| presents the minimisation procedure for the 
QMI in the general dimensional case using a trick that 
involves enlarging the orbit to its convex hull. This shift 
of perspective seems to be required in order to prove that 
the minimum QMI state is classically correlated, and no 
simpler argument has been found. However the QMI 
is not constant on this set of classical states, and us- 
ing majorization techniques we show that the minimum 
must lie in a subset of these states, most easily described 
in "Young tableau" form. In section VI F| we identify a 
partial ordering on the representative Young tableaux in 
terms of the marginal entropies, and discuss a curious 
geometrical property we call the "see-saw effect" which 
expresses the behaviour of the individual marginal en- 
tropies when transforming from one tableau to another, 
and prohibits us from knowing which tableau is optimal 
without a deeper analysis, such as that in Ref. [5]. Sec- 



tion VI G discusses special cases of classical states for 
which we can unambiguously say how the QMI varies. 



II. SETTING THE SCENE 

Let us take a moment to define some of the mathemat- 
ical objects that will appear throughout. The global or 
joint quantum state p describing a bipartite system AB 
is a positive operator on the Hilbert space "H = T-La^T^b^ 
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where dini(?^) = d, with unit trace and spectrum A = 
{^i}i=i- It has reduced states pA and ps, the states of 
subsystems A and B obtained by taking the partial traces 
of p, which are operators on Hilbert spaces Ha and Hb- 
The dimensions of the subspaces are dim(HA) — dA and 
dim(HB) = dB so that d = dAds- 

In connection with the physical scenarios we have de- 
scribed in Ref. [2], we use the quantum mutual informa- 
tion (QMI) of the state p as a measure of total correla- 
tions between its two subsystems A and B: 

I{p)^S{pa) + S{pb)~S{p). (1) 

S{p) = — Trplogp = — ^^^^AilogAi is the von Neu- 
mann entropy of state p and the logarithm is taken to 
base 2. Often we will write S{p) = H{A), H being the 
Shannon entropy, or for qubits when there are only two 
eigenvalues A, 1 — A we use the binary entropy notation 
H{\) = -AlogA- (1- A)log(l- A). The QMI is a spe- 
cial case of the relative entropy I{p) — S{p\\pA ® Pb) 
and so it is a measure of how distinguishable the state 
p is from the uncorrelated product state of its marginals 
PA® Pb- In other words, it measures the total, quantum 
plus classical, correlations in p. The relative entropy is 
non-negative S'(p||ct) > with equality if and only if 
p — a. When cr = pa® Pb this leads to the subadditiv- 
ity of entropy S{pa) + S{pb) > S{p), which ensures that 
the QMI is non-negative. The QMI is zero if and only 
if p = Pa® PB: and it is invariant under a local unitary 
transformation I{p) = I{UA®UBpUj^®Ug) or a swap of 
the states of A and B. These symmetry transformations 
define an equivalence class of states with the same QMI. 
In the following, when we, for instance, search for a state 
that gives the lowest QMI, we quote as the solution just 
one state but it is understood that the full solution is 
modulo local unitaries and swaps. 

When talking about closed evolution, which is essen- 
tially a unitary transformation on a system, the concept 
of a unitary orbit proves useful. Formally this is the set 
i 

Op = {UpU^ : VC/eSU(d)} 

and it is called the unitary orbit of p. Unitary evolution 
preserves the spectrum of a state, hence spec(t/pC/^) = A. 
It follows that the entropy (being a function only of the 
eigenvalues) of all the states in a unitary orbit is S{p). 



III. EXTREMAL CORRELATIONS ATTAINED 
UNDER UNITARY TRANSFORMATIONS 

Our motivating examples in Ref. |2| lead to the 
question of how much a unitary operation can correlate 
or decorrelate a bipartite state. In other words, we look 
for PmimPmax G Op that givc the largest and smallest 
values of I{p). 



Since S{p,nax) = S{pmin) = H{A), we can write 

HPrnax) - I{Pmtn) = ^S{pa) + ^S{pb) ■= A/ 

where AS{p^) S{p,nax,^L)-S{p^iri,t,), P- = {A,B} and 
Pm,n is a reduced state of pm, rn — {min, max}. 

We can not use convex optimisation to solve this 
because the domain of states is not convex, instead 
we search for Pmin/max along the unitary orbit sep- 
arately (later we introduce a trick that allows us to 
minimise/maximise using convexity properties). This 
means determining the reduced states that extremise 
S{pm.A) + S{pm,B), 'm = {min, max} whilst keeping the 
spectrum of the composite system fixed. 

The entropy function, when it is defined over the spec- 
tra of all states in a d-dimensional Hilbert space, is a 
concave function over the set of all possible eigenvalues 
which, being essentially probability vectors, form a con- 
vex set. It ranges from for pure states, whose eigenval- 
ues are (1, 0, 0, . . . , 0) up to a permutation, to log d for 
maximally mixed states, whose eigenvalues are uniform 
and equal to 

A. The maximum QMI state pmax 

The derivation for the maximally correlated state on 
the unitary orbit Pmax appears in Ref. [2] and will not 
be repeated here. The result for dA = dB is 

N 

Pmax =^M'^t){'^i\, (2) 

1=1 

where {I'&i)} is any generalised Bell state basis N = 
d\ and I{praax) = 2 log " H {K) . [22] 

B. The minimum QMI state pmin 

Finding pmin on the unitary orbit is more difficult than 
finding Pmax in general. One cannot fully decorrelate 
an arbitrary state via a unitary transformation hence 
I{pmin) > and, unlike the maximum QMI case, the 
marginal entropy sum S{p.rmn,A) + S{pmin,B) depends 
on the spectrum A. 

The portion of the QMI that varies on a unitary orbit 
is S{pa) + SipB) = H{Aa) + H{Ab), Af, = spec(p^), 
hence the domain of interest is the spectra Aa,Ab of the 
reduced states of all the states in Op. Calculating this 
set of compatible reduced states, given the spectrum of 
the global state, is rather formidable and is dubbed the 
"quantum marginal problem" [SHE]- 

C. The primitive case of two qubits 

The simplest possible case is when the subsystems are 
both qubits, = (is = 2, and we treat this problem 
separately. 



As a physical system it provides the smallest imagin- 
able bipartite scenario, and of course any physical imple- 
mentation of such a system would require the challeng- 
ing task of isolating it from the environment. Due to the 
small size of the systems one would expect proportion- 
ally small fluctuations to take place even for moderate 
temperatures, however as an illustrative example of a bi- 
partite quantum system the two-qubit scenario obeys the 
same laws as the general (Ia x case, but has the ped- 
agogical advantage of admitting an exact solution. 

Let us adopt the convention that the eigenvalues 
in spectrum A of the joint state p are arranged in non- 
increasing order Ai > A2 > A3 > A4 and denote the two 
eigenvalues of the reduced state as A^, 1 — A^ where 
< A^ < e {A,B}. Then pa,Pb are valid reduced 
states of a state in Op if and only if their eigenvalues 
satisfy the following inequalities [7] 

Aa > A3 + A4 (3) 

Ai3 > A3 + A4 (4) 

AA + Ai3 > A2 + A3 + 2A4 (5) 

|Aa - Asl < min{Ai - A3, A2 - A4} (6) 

A sketch of the region 7^(A) (or just TZ when there is 
no ambiguity) that these inequalities define is presented 
in figure [Tj The lines fp,gp,hp, p = 1,2, represent the 
equality conditions for these inequalities (|3| - ([6]). Fig- 
ure [2] is provided to see the different shapes that TZ can 
have depending on the rank of A (actually the qualitative 
difference in the shape of TZ depends on the degeneracy 
of A). In a slight abuse of notation, we use TZ for the 
geometric region of compatible reduced states in A^i, As 
space, the set of coordinates contained within that region 
and the states pajPb they correspond to. 

Minimising the QMI reduces to finding the minimum 
of the classical binary entropy sum H{Xa) + H{Xb) over 
TZ. It is not quite as simple as it may seem because 
the domain is not convex, see section [VI A[ We first vary 
each of the terms in the sum separately. Concavity of the 
entropy (and symmetry about its maximum at A^ = |) 
implies that H{Xfj_) decreases as A^ tends to zero. The 
line that bounds TZ closest to zero is ^pt A^ + A^ = 
A2 + A3 -I- 2A4 so the minimum of H{Xa) + H{Xb) occurs 
somewhere on this boundary of TZ. We show in Appendix 
A that the minimum of this function occurs at points 
Ri , i?2 in figure [T] where 



Xa — A3 + A4 
As = A2 + A4 



(7) 

(8) 




FIG. 1: Typical set TZ of unitary orbit marginal eigenvalues 
Aa,A_b from inequalities (|3| - with equalities equal to 
the boundary lines labelled fp, gp, hp, p — 1,2. g\ refers to 
the case when A2 = A3, light plus dark grey area, and g2 to 
A2 > A3, just the dark area (note that the lines hp, p — 1,2 
could be different for these two cases but here they are drawn 
as the same for simplicity). The marginal eigenvalues for the 
minimally correlated state pmin are situated at A when gp = 
gi and at R\,R2 when gp ~ g2. The maximally correlated 
state Pmax has marginals at B. 
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or vice versa, and the joint state (modulo the equivalence 
class of the QMI) which corresponds to it is 



FIG. 2: Regions TZ of allowed A^ (y-axis), \b (a;-axis) when 
the joint state of two qubits, with spectrum A — {Ai}*^i, 
has various ranks: A = (a) {1,0,0,0}, (b) {0.8,0.2,0,0}, (c) 
{0.5,0.5,0,0}, (d) {0.6,0.3,0.1,0} (for the entire shaded region). 
Aa,A_b G [0, |]. For each spectrum, the filled circles corre- 
spond to Pmin, the hollow ones to pmax- In (d), a state with 
energy E defines the set TZe oi states which could have en- 
ergy E. It is bounded from "above" by the thick black line, on 
which the state itself is situated. The maximally correlated 
state in TZe is at q. 



E 



X,k\e,){e,\®\Jk){fk\ 



where Aj^ is a reindexing of A^ (Aqo = Ai, Aoi = A2 
and {|ej)}, {|/fc)} are qubit basis states for systems A and 



(9) 



) 



B respectively. The QMI for two qubits has a minimum 
value of 

/(p™„) = H{X^ + A4) + H{X2 + A4) - H{K). 
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Here we have used the binary entropy notation: H{X3 + 
A4) = H{\i + A2) = H{\i + A2, A3 + A4). 

Any state p — Ai|V'i)(V'i| diagonal in the basis 

\4>i) along the unitary orbit of Pmin can be transformed to 
Pmin by the unitary Umin = Y^t=i where \efi) = 

\ej)\fk) are products of qubit basis states and j,k = 0, 1 
are the binary representation of i. 

So to answer the initial question, for a two-qubit sys- 
tem, the maximum that the QMI can change by over 
the domain of states in Op with spec(p) = {Xi}j^i, 
Ai > Ai_|_i, is 

A/ = 2-iJ(Ai+A2)-i?(Ai + A3). (10) 



H{Xa) + H[E — A^), and its maximum occurs when 
Xa — Xb = ^T point q in the figure. Hence the max- 
imum possible initial value of the QMI for a state with 
total average E is 2H{^). 

After the global unitary transformation UpU^ = 
p' , the local subsystem states evolve to = (1 ^ 
X'^)\v)^{v\p + X'^\v)^{v\^ with a; < i and {\v)^,\v)^} 
is a qubit basis 

|i;)^ = cos^|0)-e-^^-sin^|l) (13) 
|i;)^ = sin^|0)+e^^''cos^|l) (14) 



D. The Imposition of a Global Constraint: Energy 

Conservation 

In this section we seek the maximum change in the 
QMI 

AIe - I{p') - lip) = S{p'a) + S{p'b) - S{pa) - S{pb) 

(11) 

subject to p' = UpW where U is now an energy- 
conserving unitary satisfying E = Ti^pH) = Ti{p'H). 
H is a Hamiltonian for the composite system AB, taken 
to be local H = Ha + Hb , and the initial reduced states 
are thermal with respect to Ha , Hb ■ These are the as- 
sumptions for the heat flow model between two quantum 
systems. Example 2 in Ref. Below we find the maxi- 
mum AIe when A and B are qubits. 

We take the local Hamiltonians Ha,Hb of the initial 
subsystem states to be equal [53] : the ground and excited 
states energies are zero and one respectively and the en- 
ergy eigenbasis is the computational basis {|0), |1)} thus 
Hp = \l){l\p. We can write the initial (thermal) reduced 
state of p as 

p, = (i-a;)|o)(o| + a;|i)(i|, 

where A^ — e^^/-^f'/Tr(/9^) < 1 and Tp is the tempera- 
ture of system p. 

Let us suppose that during the energy exchange inter- 
action (a global energy conserving unitary on p) correla- 
tions decrease, i.e. the two systems A, B, though locally 
thermal, are initially correlated. We would therefore like 
to minimise S{p'a) + S{p'b) subject to total average en- 
ergy conservation. The initial total energy set by the 
initial state p is 

E = Tr{pH) = Tt{paHa) + MpbHb) = A^ + A^ (12) 

with < £; < 1. 

Note: The point (A^, A^) G 7?. is on the line A^ + A^ = 

E, the thick black line in figure [2] (d). The joint states 
represented by this line all have total energy E and 
their reduced states are thermal. The QMI, or equiv- 
alently just H{Xa) + H{Xb), is concave on this line, 
which is easily seen when the function is written as 



The final total average energy E' — Tr{p'^HA) + 
Tt{p'bHb) [patPb are not constrained to being diago- 
nal in Ha^ Hb) is 

E' = 1-^(cos6'a + cos6'b) + A^cos6'a + AbCos6'b (15) 



In order for energy to be conserved, E' — E. This con- 
dition and equation (15 1 gives a collection of straight 



lines, defining a region TZe Q Ti- representing joint states 
which could have total energy E for some choice of angles 
9a, 9b- 

TZe given by A^i + Xb < E, for the derivation see 
Appendix B and for a visualisation see figured (d). The 
energy conserving region always contains the minimally 
correlated state, hence the greatest possible change in the 
QMI at constant average energy E occurs when the initial 
state occurs when A^ = A^ = ^ (these marginals could 
define multiple global states, even after factoring out the 
states in the equivalence class of the QMI, see section 
IV A) and the final state is the minimally correlated one, 
Pmin given by equation ([9|, with 9a, 9b chosen such that 



equation ( 15 1 is satisfied (there is a range of angles which 
do this). We have then that 



A/b = 2i/ ( I ) - (Ai + A2) - if(Ai + A3) (16) 



The way to think about the region TZe intuitively is 
that it contains final states whose marginals in general 
have longer Bloch vectors than the initial p, achieving 
this higher local purity is made possible by using up the 
correlations in the joint state. The constant energy con- 
dition is met by rotating the Bloch vectors off the energy 
axis by the angles 9 a, 9b, thereby adding coherence what 
was originally a thermal state. 

Note that if we insist the also the final states are ther- 
mal, then the energy conserving region is just the set of 
states in TZ that are on the line Xa + Xb = E. The max- 
imally correlated state is then still at A^ = A^ = ^ but 
the minimally correlated one is at the edge of the line. 
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IV. ANALYSIS OF THE SET OF REDUCED 
STATE SPECTRA 7^ 

We appealed to the collection TZ of reduced states com- 
patible with a joint state having a fix spectrum A in order 
to find the maximally and minimally correlated states, 
with and without an energy conservation constraint. It 
would be useful to know how the action of a global uni- 
tary acting on the joint state translates to dynamics in 
TZ, but first we must tackle an issue of degeneracy. 

A. One point refers to a range of compatible joint 
states 

A two qubit density matrix is specified by 15 parame- 
ters, a point in TZ defines only 5 independent parameters 



(Aa, As, Ai, A2, A3), hence there is a lot of degeneracy as 
each point in TZ corresponds to a whole set of compati- 
ble joint states. Modulo the parameters that don't affect 
the QMI (the directions of the Bloch vectors, which re- 
quires 4 more parameters) there are still 6 parameters 
that can be freely chosen. In appendix A we show that 
the coordinate in TZ that represents p„iin is given by a 
unique state (modulo local unitaries and swaps). How- 
ever there are points in TZ for which this is not true, these 
points determine states that have the same QMI but are 
not related by local unitaries or swaps, hence they would 
have different amounts of correlation, or entanglement, 
if we used another measure that, say, depended on the 
correlation- or T-matrix of the state, whose elements are 
tij = Tr((T' (g) (T^/o), i,j = 1, 2, 3 and cr* is a Pauli matrix. 
We illustrate this with an example of a "triple point" in 
7^(A). Defining 



p(a, P, 7, S, cos 9, cos 0) := 



a + 5 + {a — S) coscj) 
11 
2 I 

{S ~ a) sin (p 





+ 7 + (/? - 7) cos ( 
(7-/3) sine 




((5 — a) sin 

(7-/3) sine 
/3 + 7- (/3-7)cos6l 

a + 5 — {a — 5) cos cj) 

(17) 



where the matrix is written in the computational basis, 
three states that correspond to the same Xa,^b but 
are not connected by a symmetry of the QMI are 
p(Ai, A2, A3, A4, ^^), p(A2, Ai, A3, A4, 0, ^ 



-A4 

^^4' A1-A4 

p(A3,Ai,A2,A4,0, 1) and 



Xa = Ai 



-(1-(A3-A4)). 



(18) 



We know also of "double points" and correspondingly 
pmin is a "single point" , shown in Appendix A. 

Despite the value of the QMI not being sensitive to 
these different states, their dynamics under energy con- 
serving unitaries leads to different values of A/, and it 
is precisely the change in QMI which is the physically 
relevant quantity, see Ref. [2J : it is in this way that these 
hidden correlations are revealed. 



B. Action of unitaries on p 

There are many ways in which we could analyse the ef- 
fect of a global unitary on a two qubit state, with regards 
to its position in TZ. To simplify the task, we concentrate 
on some "elementary" transformations only rotating in 
two-dimensional subspaces and discuss what their conse- 
quences are for energy conservation. 

Let us assume that the local bases are fixed in the com- 
putational basis, then the individual and total Hamilto- 
nians are Ha = Hb = |1)(1|, i? = Ha®1 + 1® He ^ 



|01) (01| + |10)(10| +2|11) (11|. There are two key unitary 
transformations (matrices are all written in the computa- 
tional basis) that are effective two-dimensional rotations: 



exp( 

/I 




\ 

cos(§) sin(|) 

-sin(f) cos(|) 

1/ 

exp(-iij('=V/2) 

/ cos(|) sin(|) \ 

10 

10 

V -sin(f) cos(|) / 



U^°\e) = exp(-iij(°)6l/2) (19) 

(20) 

\0 1/ 

C/('=)(0) = exp(-iij('=V/2) (21) 

(22) 

where the interaction Hamiltonians are H^°^ = ^{X ® 
Y - Y ® X) and H^^'^ = -\{X ®Y + Y ® X), X,Y 

are Pauli matrices. Since is a strong 

energy conserving unitary (SECU). In terms of its action 
on the reduced state eigenvalues (and therefore on the 
QMI when starting from states that are locally diagonal 
in Haj Hb which is all that we require, see example 2 in 
Ref. [2]) it is the most general unitary that commutes 
with the total Hamiltonian. 

The matrix U^°\e) rotates in the energy degenerate 
"odd" subspace {|01),|10)} and U'^''\(j)) rotates in the 
"even" subspace {|00),|11)}, but is not energy conserv- 
ing. U^°\9) and C/'-'^' (0) commute. Let the initial state 
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be diagonal in the computational basis 



Po = 




(23) 



and the elements on the diagonal are some permutation 
of eigenvalues A = To see how the unitaries 

transform po in the region TZ we need to look at the 
(lowest) eigenvalue of the reduced states. A state evolved 
under these unitaries looks like p(q;, /3, 7, i5, cos 0, cos 0) 
in equation (17 1, it is often called an "X-state". Note 
that p{a, (3, 7, 6, cos 6, cos (f)) has reduced states that are 
diagonal in the local bases Ha, Hb, that is, its local states 
are thermal (the spectrum can be equated to a Gibbs 
distribution). In this notation po = picn, /3, 7, S, 1, 1). 



The reduced state 
p{a, j3, 7, 5, cos 9, cos (j)) are 



eigenvalues 



AA(a,/3,7, (5, cos6',cos0) = 
(1 — (/? — 7) cos 9 — {a — 5) cos 4>) 

Xsioi, f3,j,S,cos9,cos(j)) = 
- (1 + (/? - 7) cos 9- (a-S) cos 0) 



of 
(24) 

(25) 



if we pick min{a — (5, /3 — 7} — /3 — 7. Care must be 
taken to ensure that these eigenvalues are the smallest 
depending on the ordering used for a, /3, 7, S. 

Following from this, Xa + = 1 — (a — (5) cos <j), hence 
varying (f) provides a set of lines of the form A^+A^ — Ci. 
On the other hand. As — A^ = (/? — 7) cos 9, thus 9 defines 
the lines Xa — Xb + C2. Therefore, when the initial 
state is that given in equation ( |23[ ), we can describe the 
evolution of the odd and even unitaries fairly simply in 
terms of the marginals and the regions that they trace 
out in 7?,. It is easy also to see their effect on the QMI 
(in simplified notation): 

/(cos 9^x)^ H{c - mx) + H{c + mx) - H{A) (26) 
/(cos (/) = y) = H{a - ny) + H{b - ny) - H{A) (27) 

where a, b, c, m, n are all constants related to terms in 



equations ( 24 ) and ( 25 ) , and the QMI is (separately) con- 
cave over the domains x and y. 



It is shown in section VI D that of the 4! possible or- 
derings of Xi, i = 1,2,3,4, there are only three classes 
of permutations which give distinct values of the QMI. 
Representative states for each class are taken to be 
tW = p(Ai,A2,A3,A4,l,l), r(2) = p(A2, Ai, A3, A4, 1, 1), 
T^'^) = p(A3, Ai, A2, A4, 1, 1). Note again that, especially 
when calculating energies, care must be taken so that 
the reduced states of each of these are in the form 
P,. = (1 - A^)|0)(0| + A^|l)(l| with < A^ < i - lo- 
cal unitaries may be used to ensure this condition. Using 
each of these as an initial state pq, figure |3] depicts the 
effect of U^'^\ U'^°^ on the each of their marginals. 
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FIG. 3: (Color online) The reduced state eigenvalues of 

U'^°\e)poU'-°\e)\ ^7('='((^)po^7''=^0)^ O < 6»,0 < tt when 
po = r'"', i = 1, in blue, i = 2, in red and i = 3 in green, 
r''^ is a classically correlated state whose marginals are also 
plotted and a triple point TP is shown when three states with 
different T-matrices but same values of QMI coincide. The 
global spectrum is Ai = 0.6, A2 = 0.3, A3 = 0.1, A4 = 0. 



We can clearly see the triple point described in section 
|IV A[ labelled TP in figure [3] as well as many other over- 
laps corresponding to triple and double points. A corol- 
lary of this is that starting at a triple/double point, the 
unitary f/(°) generates portions of the line A^ + A^ = Ci 
that are of varying length, depending on which exact 
state we started with. The longer the length, the lower 
the correlations can go since the QMI is concave on the 
line. 

Comparing figure [3] with the full allowed region 7?.(A), 
see for instance figure [2] (d), we see that just these two 
unitaries can evolve the classically correlated state to al- 
most every point in TZ{h). The "wings" that are missing 
in figure [3] can be obtained using one more unitary 



u{0 



exp(-i//e/2) 

/ cos(|) sin(|) 
— sin(|) cos(|) 



v 



q\ 



1 

1/ 



(28) 



with H ^ ~^(I<SiY + Z®Y). This unitary acting on 
p(Ai, A2, A3, A4, cos^, 1) results in figure [ij By consid- 
ering the transformation this unitary generates on the 
marginal eigenvalues, it is simple to show that these three 
unitaries generate all points in TZ. To summarise, we 
have identified how to reach all points in the allowed re- 
gion TZ starting from classically correlated states using 
just three unitaries, overlaps are easily seen to be dou- 
ble/triple points and the joint state they correspond to 
is known. 

Let us take a closer look at unitary evolution for states 
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FIG. 4: (Color online) Same as figure [s] with the marginal 
eigenvalues of C^(^)r'^^f7(^)^ added, < ^ < tt, this now is 
equal to 7?.. 



in the energy conserving region TZe C TZ, defined as the 
set fulfilling Aa + As < E. We set the initial state now 
as the maximally correlated state for this subregion po = 
it is at the point on the Une A/i + As = E where 
Xb- The local and total Hamiltonians are still the 



Pe 



Xa 

same (diagonal in the computational basis) . To conserve 
energy, the unitary evolution of po must obey 

E = TiipoH) = TiiUpoU^H) 

There are two types of unitary that satisfy this: SECUs 
Us that commute with the total Hamiltonian (above we 
stated that the only SECU for this Hamiltonian is [/'°') 
and weak energy conserving unitaries (WECUs) Uw that 
do not commute with the total Hamiltonian but preserve 
the total energy on average. 

Depending on the exact form of p^"*^ (could be e.g. a 
triple point) the SECU evolves it along some portion of 
the line Xa + Xb = E. Getting to other states in the 
region TZej for instance where Xa + Xb < E, and con- 
forming with the energy conserving condition, equation 
(151, requires WECUs, which can shrink Xa,Xb, and lo- 



cal unitaries to rotate the local bases by appropriate an- 
gles 9a, Ob- This generally amounts to elongating the 
local Bloch vectors and rotating them off the energy axis 
in order to preserve total average energy. Consequently, 
the reduced states are permitted to be non-diagonal in 
the energy eigenbasis and so they are no longer thermal, 
however we argue in Ref. and the next section that 
this is still physically and thermodynamically relevant. 
Note there is a range of angles 6a , Ob that satisfy equa- 
tion ( |15[ ) and the amount energy exchanged between the 
qubits depends on these angles. 

There is an interesting conclusion to be noted from all 
of this: sometimes a SECU can only transform a state 
with total energy E along a portion of the line Xa + Xb = 
E but there are states that are on Ayi -I- Ab = E mTZ that 



are not reachable by the SECU and these unreachable 
states are locally thermal. This means that a WECU is 
required to get to these states. Hence if we, say, restricted 
our heat flow model (in the next section) to initial and 
final states being thermal, with total energy conserved, 
it would not be general enough to only consider SECUs, 
one would still have to include the WECUs for the full 
analysis. 

Another property of transformations under unitaries 
that commute with the Hamiltonian is that the quantum 
variance S — Ti[pH'^] — (Tr[pi/])^ remains constant as 
p — )■ Us{t)pUs{ty ■ We assume constant average energy 
Ti[[pH] = E so changes in S only depend on the first 
term. All moments of H, Tr[pH''] for k = 1,2... are con- 
stant for SECUs. This means that the statistics of all 
the measurements of the energy of Us{t)pUs{t)^ will be 
the same, due to total constant energy and the reduced 
states being diagonal in the energy eigenbasis. Quantum 
uncertainty due to coherence is introduced when the ob- 
servable and the state are not locally codiagonal, and in 
general the term Tr[pi?^] will increase for unitaries the 
do not commute with H . 



V. APPLICATIONS TO PHYSICAL PROCESSES 
WITHIN A CLOSED QUANTUM SYSTEM 

Having developed a body of results in the previous 
section we apply this work to several scenarios that arise 
in quantum thermodynamics [5], such as equilibration 
and heat flow between two initially thermal states. As 
before, the underlying principles and equations are quite 
general, and for simplicity we make extensive use of the 
illustrative two-qubit case, and highlight points where it 
differs qualitatively from the more general case. 



A. Heat Flow 

The model for the heat flow, example 2 in Ref. [2], 
between two systems A and B is based on three assump- 
tions: (1) the global evolution of the composite system p 
is unitary p ^ p' = UpW , (2) the total average energy of 
p is conserved, E = Tr(pi?) = Tr(p'7?) [M] and (3) the 
reduced states pa, Pb of the initial state p are thermal. 

This model places a limit on the amount of heat that 
can flow from one system to another. It adheres to con- 
ditions (l)-(3) above and, using the fact that the free 
energy is minimised by a thermal state, yields the in- 
equality [2j 



1 



kTA 



kTr, 



> Ah 



heat ; 



(29) 



where Qa is the heat into system A, Ta (Tb) is the ini- 
tial temperature of system A (B), k is the Boltzmann 
constant and A/jjoat is the difference in the mutual infor- 
mation before and after the heat exchange interaction. 
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its greatest value for two qubits is displayed in equation 

PI, |A/hcat| < A/b. 



B. A Digression on Physical Assumptions 

To reduce misconceptions, we would like to discuss 
several key points about the physical assumptions un- 
derlying the application of our formal results to physical 
scenarios involving closed quantum systems. We use the 
heat flow model as a base for our discussion, but bear in 
mind that it is not limited to this one situation, nor is it 
restricted to qubit systems. 



Range of Applicability 



The first issue that one may raise about inequality ( 29 ) 



is on the usage of the free energy and the definition of 
temperature. Only the initial temperatures of the sys- 
tems A, B appear in (29), and the temperatures of the 



final states may be different or may not even be defined 
(they are not restricted to be thermal states). If this is 
so, how can the free energy be used here? The answer is 
that we employ the free energy functional Fh.t{p) which 
is initially an information theoretic quantity and, as a 
mathematical object defined over all state space, can be 
applied to any quantum state, not just Gibbs states. The 
Gibbs, or thermal, states p*'* — exp{—H/kT)/Z are spe- 
cial in as much as they are the states that minimise the 
free energy functional 

FH,T{p)^T^4Hp]-TSip), 

S{p) is the von Neumann entropy of p. Fh,t{p*^) co- 
incides with the usual thermodynamic free energy re- 
sult. In fact a change in the free energy functional 
Fh,t{p) ~ Fh,t{p^'^) is proportional to the relative en- 
tropy S{p\\p^^), which has an operational meaning in 
terms of the distinguishability of two states, and empha- 
sises the physical relevance of the free energy functional. 
The positivity of the relative entropy for all quantum 
states in turn corresponds to a directionality for any ther- 
modynamic process: AFh.t > 0. In standard situations 
this corresponds to the free energy being minimised by 
an equilibrium state. The heat flow goal is then to ex- 
plore the asymmetry in hot — ?> cold and cold — ?> hot heat 
flow from Eq. ( 29 ) , which is a broad technical constraint 



dependent only on initial temperatures and does not rely 
on defining final temperatures. 



2. Local Thermality and Entanglement 

The free energy functional constrains any changes in 
local variables, irrespective of coupling strength, and it 
is the physically relevant quantity that provides the di- 
rectionality for such internal statistical processes. Local 



concepts, such as the initial temperatures of A and B, 
can still limit non-local properties: the reduced states 
are initially entirely indistinguishable from thermal states 
(when restricted to doing only local operations on them) 
but the presence of correlations means that we do not 
assume "molecular chaos" [9]. When one works in this 
broader framework we still have directionality for any 
energy conserving process which comes from the posi- 
tivity of the relative entropy. The dependence on tem- 
peratures just reflects the assumption of initial thermal 
marginals but also places a limit on the amount of cor- 
relations possible (for example two systems which are lo- 
cally "too cool" cannot be maximally entangled), which 
in turn bounds any anomalous heat flow. 

It is very much permissible to have a bipartite system 
whose marginals are thermal but the global system is al- 
lowed to be correlated or even entangled. Indeed recently 
it has been shown [TUl that local thermality arises typi- 
cally for constant energy pure states. The degree of corre- 
lation is also something we, a priori, do not limit, however 
any correlation in p must be consistent with its reduced 
states PA, Pb being (initially) Gibbs states. This con- 
dition means that outcomes of local operations on, say. 
Pa must be indistinguishable from doing the same oper- 
ations on an uncorrelated thermal state with the same 
temperature and Hamiltonian. The correlations amount 
to extending the description of the system from local ob- 
servables {Ha, Hb} to {Ha, Hb,Ha ® Hb}- 



3. Bipartite Interactions and Heat Flow 

Let us clarify what it means when we say systems A 
and B interact. Interacting here means that the systems 
undergo a global unitary evolution, which conserves the 
total energy. This evolution could be thought of in dif- 
ferent ways. It could be controlled by an experimenter 
who switches on an interaction V(t) for some finite time. 
Or we can think of it as a scattering process: initially 
the particles are assumed to not interact, e.g. they are 
far apart, hence Hi = Ha + Hb, but there could al- 
ready exist an unknown degree of correlation between 
them. The role of the experimenter would then be to 
bring the particles together so the particles can interact 
via some natural entropy and energy conserving process. 
The only action the experimenter takes is in bringing the 
particles together and then apart, controlling the inter- 
action strength (how close they are) and the interaction 
time. Both these circumstances are modelled as: To- 
tal Hamiltonian H = Ha + Hb for t < 0, t > T, and 
H = Ha + Hb + V{t) for < t < T, and no net work 
is done on the total system by assumption. Note there 
is no explicit time t in our model but we assume there 
is some "clock" that governs the temporal length of the 
interaction. 

During part of the total interaction time T, the ex- 
perimenter may have to put in or extract energy out 
of the system in order to controUably interact A and 
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B, however the total system AB is required to finish 
in a state such its total energy is equal to its initial 
total energy, and so overall the experimenter does no 
work. In this way we can formulate energy conserva- 
tion in our heat exchange model as a zero work process: 
W = TT[Hpfj,g\-TY[Hp\g\ =0,H = Ha+Hb. Energies 
exchanged between the subsystems are assumed inacces- 
sible for external work. As a result we call the internal 
exchanged energy "heat" . This is in turn supported by 
the fact that the local entropies change. This is still a 
valid measure of average energy even when, say for sys- 
tem A, [p;^, Ha\ 7^ 0, because energies are measured with 
respect to the original energy eigenbasis defined by Ha 
even when the final states are no longer diagonal in this 
basis. This is validly done elsewhere in the literature, see 
for instance Refs. [TT] and [T^]. 

Normally heat flowing from cold to hot is associated 
with an external source of energy or heat pump driving 
the process, however this can also occur within a closed 
system by using up internal correlations. In this sense, 
correlations are a resource. The model is in a way an 
embodiment of the Clausius form of the second law: "No 
process is possible whose sole result is the transfer of heat 
from a body of lower temperature to a body of higher 
temperature " . 

One usually considers thermodynamic processes from 
an "outside view" where heat baths and work reservoirs 
externally control the system. However one can switch 
to an "inside view", where external influences are de- 
scribed as eff'ective changes in the system's Hamiltonian. 
This leads to an alternative "inside view" definition of 
heat and work [13j . The unconventionality in our heat 
exchange scenario arises from the fact that we do not as- 
sume interactions with heat baths or work reservoirs, all 
of our interest lies in the internal view in order to ob- 
tain as general a set of constraints on the evolution as 
possible - i.e with as little model dependence as possi- 
ble. We restrict ourselves to a local thermodynamic level 
of description, the presence of correlations manifests it- 
self in negative heat flow, inferred from changes in (Ha), 
{Hb}- The actual unitary that governs the interaction of 
the systems A and B may be unknown to the observer. 
Thus the internal view in |13| could in principle be ap- 
plied here, but this is not our regime of interest, since 
the observer would have to know the initial correlations 
between A and B and the type of unitary used. 

Stationary heat flow is normally maintained if each sys- 
tem is coupled to a large external heat bath leading to 
dissipation and irreversibility. If this were the case then 
the effect of correlations might be washed out. How- 
ever we work in a different regime: we do not want to 
assume an open system, there is no environment, no dis- 
sipation and we focus on internal processes of closed sys- 
tems. Energy exchange occurs only between A and B 
(which in principle can be any size) and it is controlled 
by the interaction "time". Once A and B stop interact- 
ing, they are static and energy exchange desists. It is in 
this way that we obtain heat exchange - not by dissipa- 



tion. Since A and B have relaxed after the interaction, 
they can be assumed to be "in equilibrium" (but not nec- 
essarily Gibbs states) since their states may (see below) 
remain stationary. 



4- Strong versus Weak Coupling 

We now turn to the more subtle issue of coupling 
strength. It has been shown that without weak coupling, 
serious problems arise [HI [T3] such as apparent local vio- 
lations of Landauer's principle. We make no assumptions 
about coupling strength in our model, nor do we need to. 
First we must determine what exactly is meant by cou- 
pling strength in our situation; in interacting systems 
there can be a notion of internal and external coupling 
strengths. Usually coupling strength is a term associated 
with open systems, i.e. the interaction of a system with a 
thermal bath, the total Hamiltonian of the system being 
H = Ha + Hb + gHint (t) , so we can make the identifi- 
cation V{t) = gHint{t), where g is the coupling strength. 
Roughly speaking, in an open system approach, A might 
be low dimensional and B a bath (infinite dimensional). 
When g is small (weak coupling) the subsystem of inter- 
est, A, is assumed to stay diagonal in its original local 
energy eigenbasis, but in the long time limit thermalises 
to a Gibbs state having the temperature of the bath. 
When g is large (strong coupling) the joint system AB 
evolves to a Gibbs state, with Hamiltonian H, and A is 
no longer diagonal in Ha- 

As mentioned, we do not have an interaction "always 
on" but consider a transitory interaction (scattering). 
Before and after, H = Ha + Hb and no ambiguity arises 
in the specification of subsystems and their individual 
properties. As such we need not impose restrictions on 
the transitory V{t). A useful comparison might be in 
the formulation of thermodynamics by Lieb & Yngva- 
son [16j . In formulating the principle of entropy increase 
from the second law, Lieb & Yngvason define a thermo- 
dynamic process via "adiabatic accessibility" . The ther- 
modynamic entropy is defined as (amongst other things) 
an additive function that is non-decreasing during these 
transitions. They require that the initial and final states 
of the system be in equilibrium, in the sense that the 
system has relaxed to a stationary state but they make 
it very clear that the process in between need not be 
slow nor gentle, in fact it can arbitrarily violent (e.g. see 
Lieb & Yngvason's "adiabatic gorilla"). Hence there is in 
principle no restriction on the coupling strengths under 
these asymptotic conditions. 

One might however worry about internal coupling 
strengths in the heat flow model. There, weak coupling 
would imply that the final states after the unitary inter- 
action are still approximately locally diagonal in Ha, Hb 
whereas strong coupling corresponds to final states that 
are not diagonal in Ha, Hb- We do not assume, a priori, 
any particular coupling strength; the governing equation 
( 29 1 is independent of coupling strength. The systems 
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A and B can potentially go to thermal or non-thermal 
states after the interaction because, as mentioned above, 
the free energy functional that governs the directional- 
ity of the process is defined over all state space (not just 
Gibbs states) , what matters is only that the initial states 
are locally thermal. 

In this sense such a coupling concern is in allowing 
the final reduced states to not be Gibbsian, leading to 
[P'^j^,Ha] ^ and similarly for subsystem B. However 
this is not necessarily a problem when we remember that 
the total system AB is isolated, i.e. after the interac- 
tion, there is no thermal bath for it to equilibrate with 
and once the particles are separated after scattering, the 
reduced states have fixed energies with perhaps coher- 
ent oscillations in their energy eigenbasis (even though it 
might be the case that we cannot ascribe any notion of 
local temperature - in fact we don't need to if interested 
solely in questions of energy exchange). What is differ- 
ent in the strong coupling case is that the von Neumann 
entropy does not equal the thermodynamic entropy any 
longer (if one defines the thermodynamic entropy to be 
for the state having energy Tr[pi7] but which maximises 
von Neumann entropy, defined thus in jllj). Despite 
this, in the event of such non-static outcomes, familiar 
entropic and thermodynamic behaviour is still observed 
with respect to the initial equilibrium states. This be- 
haviour stems once again from the directionality implied 
in Eq. ( 29 1 , itself born from the positivity of the relative 
entropy with respect to the initial Gibbs states. 



C. Towards a Generalised Collision Model 

In an attempt to explain equilibration between two sys- 
tems, Partovi [T7j devised a collision model between two 
particles. We outlined its principles in example 3 in Ref. 
[2] and claimed that when the two systems are qubits, 
one of Partovi's restrictions can be weakened. Here we 
elaborate on why this is true and furthermore we explain 
why in the higher dimensional case, Partovi's initial as- 
sumptions must hold. 

The fundamentals of the collision process is sum- 
marised here: initially the two systems are described by 
a joint product state p, they interact via a global energy 
conserving unitary and then they decorrelate, forgetting 
any correlations that are produced and rendering the pro- 
cess irreversible (this moves the system to a new unitary 
orbit). More compactly: p = pj^® pg p' — UpW -> 
p'a ® p'b- the end of these steps, the local entropies 
and global entropy has increased, H^Sa + AS's > 0, 



AS* > 0. This whole procedure is then reiterated un- 
til the entropies are maximised and the systems reach 
equilibrium. 

We propose a modified version of the model where, in- 
stead of the systems completely decorrelating after the 
unitary interaction, the joint state of the system de- 
phases to a classically correlated state by losing its "off- 
diagonal" components (again this process transfers the 
system to a new unitary orbit). Note that this requires 
the introduction of local bases. This process may be more 
physically meaningful because dephasing occurs faster 
than decorrelation in real systems. 

Let (7, a' be two minimum QMI states on some unitary 
orbits. Then does the process a —y p' = UaW — > a' 
achieve the same entropic increase as the original Partovi 
scheme? a' is the dephased p' . 



1. Two qubits, cLa ~ ds = 2 

The total Hamiltonian for the system can be taken 

as H = Ha + Hb = |1)(1U ^ Ib + U ^ |l)(l|s = 
|01) (01 1 + 1 10) (10| -I- 2| 11) (11 1 . It has an energy £; = 1 de- 
generate subspace spanned by {|01), |10)}. The most gen- 
eral form of the unitary that commutes with this Hamil- 
tonian (in a loose notation) is Uh = ^e=o+Ue=i+^e=2, 
where H^^o = |00)(00|, H^^s = |11)(11| and Ue=i is a 
general 2x2 unitary rotating in the energy degenerate 
subspace. Without loss of generality (in terms of action 
on the marginal eigenvalues), Uh need only rotate in a 
two-dimensional plane of the effective {|01), |10)| Bloch 
sphere. Thus Uh is given by C/(°) in equation [m. 

The initial minimally correlated state is cr = 
Ai|00)(00| + A2|01)(01| + A3|10)(10| -t- A4|ll)(ll|. It is 



located at the minimally correlated vertex of its reduced 
state spectra polygon TZ{a), c.f. where Pmm is located in 
figurejl]- at i?i or i?2. Its energy conserving region TiE{(y) 
therefore is just the line Xa + ^^b = E = A2-I-A3-I-2A4, i.e. 
the lines /i or /2 in figure[l] and so Uh (a SECU) evolves 
a to all states in TZe{o')- The fact that [Uh,H] — 
means not only that energy is conserved, but also that 
the reduced states will remain diagonal in their energy 
eigenbases and, since they are qubits, they will remain 
thermally distributed (it is always possible to identify 
the qubit eigenvalues with a Gibbs distribution with re- 
spect to Ha and Hb) and so one can always define local 
temperatures. 

The first step in the collision process is cr — > p' = 
UncrUlf and the final state p' ends up somewhere along 
the constant energy line forming 7?,£;((t) and looks like 



/ Ai 0_ \ 

1 pA2 + (l-p)A3 -VM1-P)(A2-A3) 

2 -v/p(l-rt(A2 - A3) (l-p)A2+pA3 

V A4 / 



(30) 
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with p = cos^ 6. Now decoherence occurs in the com- 
putational basis on the global state and the final state 



/ Ai 

pA2 + (l-p)A3 

(l-p)A2 

\ 



pAa 



(31) 



The dephasing process moves the state of the system off 
the original unitary orbit because S'(cr') > S{p'). To show 
a' is still a minimally correlated state in its orbit we must 
check the ordering of its eigenvalues. It is readily seen 
that if Ai > A2 > A3 > A4 then, since < p < 1, it is 
true that Ai > pX2 + (1 — _p)A3 > (1 —p)X2 +PA3 > A4 or 
Ai > (1 — p)X2 + pXs > pX2 + (1 — p)A3 > A4 and both 
of these orderings are valid minimally correlated states 
(one ordering is obtained from the other by a swap of 
the states of A and B which is a symmetry of the QMI) . 

The combined transformation on the state can be writ- 
ten as a stochastic map : cr — > J2t=i^kUHcrUljIlk, 
where 11^ is a projector onto the joint basis liAjs), 
k = (lAjjE) and iaiJe = 0,1. After one collision, the 
map on the initial classically correlated state a trans- 



forms it to the state a' in equation (32). After n colli- 
sions, the resulting state is 



/Ai 



V 





An) 

^2 








An) 





A4/ 



(32) 



where 



^ p(n-l)A("-l) + (1 - p("-l))A(""^) 



a(") = (1 -p("-i))A("-^) +p("-i)A("-^) 
It is simple to see that the intervals satisfy 



,(ri) ,(n) 
^2 ' -^3 



>(n-2) 
1 ) -^2 



e •■• e [Ai,A2] 



hence after each collision, the difference between 
the middle two eigenvalues Aj""* , Ag"-* decreases and 



lim„_j.oo IA2"'' — Ag"^! = 0. Also the ordering of the eigen- 
values on every new orbit (after dephasing) preserve the 
minimally correlated form. Thus after very many colli- 
sions, the local and global entropies are maximised, the 
system's reduced states have the same spectra and so 
equilibrate to the same temperature. This means that it 
is not necessary for qubits to forget the correlations gen- 
erated by the unitary interaction; if they only dephasc to 
a classically correlates state, they will still equilibrate to 
the same temperature. 



2. The qubit-qutrit case: (Ia ~ 2, cIb ~ 3 

This property does not unconditionally carry over to 
higher dimensions, which we now demonstrate by looking 
at the dA = 2,dB = 3 case. 

The Hamiltonian of the total system now is H = 
|01)(01| + |10)(10| + 2(|11)(11| + |02)(02|) + 3|12)(12|. 
Now there exists an energy E — I degenerate sub- 
space spanned by {|01), |10)} and a.n E = 2 subspace 
spanned by {|11),|02)}. Accordingly, the most gen- 
eral form of the unitary which commutes with H is 
Uh = TIe=o + Ue=i + Ue=2 + ^E=3 where IIb^o = 
|00)(00|,n£;^3 = |12)(12| and Ue=i,Ue=2 are 2 x 2 uni- 
taries rotating independently in their own energy degen- 
erate subspaces. However unlike the two qubit case, these 
energy conserving rotations are not "symmetries" of the 
QMI, in the sense that after the unitary operation and 
the dephasing, the system does not return to a mini- 
mally correlated state. For example, let the initial state 
be CT = Ai |00) (00| + A2 101) (01 1 + A3 102) (02| + A4 1 10) (10| 
A5|ll)(ll| + A6|12)(12|, with Ai > A2 > A3 > A4 > A5 > 
Ag; we pick a spectrum A = such that this is 

the its minimally correlated form, see section VI A per- 
mitted energy conserving transformation in the first step 
of the collision process is one whose only action is to 
swap A2 and A4 in a. After dephasing (which in this 
case has no effect) the state is now p' which is equal to p 
up to a permutation of the eigenvalues on the diagonal. 
This particular permutation raises the QMI of the state 
I{p') > -'^(c) (see section VIB 1 ), hence p' is no longer a 



minimally correlated state on its orbit (which here is the 
same orbit as a). This means that in the next iteration 
of the collision process one can find an energy conserving 
unitary that decreases Sa,Sb (for example, just the in- 
verse of the unitary that was initially applied to a) which 
goes against the assumptions of Partovi. 

To summarise, if we are to concur that systems tend 
to equilibrate via the Partovi collision model, where at 
every step the local entropies are non-decreasing, then 
(apart from two qubit systems) the two colliding systems 
must decorrelate to a product state after every interac- 
tion. In systems with ci > 4, swaps of energy level popu- 
lations that preserve energy do not preserve correlations. 
A corollary of this is given a Hamiltonian H = Ha + Hb 
and an energy E = Tt{pH), the set of classically corre- 
lated states {p} that are diagonal in H with energy E 
have a range of QMI values AIe,h- This is a special case 
of the energy conserving condition analysed in section 



HID above, but restricted to classical states. 



D. Global demons and local observers 

Our analysis of the correlations in a bipartite state 
can be phrased in terms of an paradoxical scenario in 
which a global Maxwell demon can confuse an observer 
Charlie who can only measure local observables. Char- 
lie is handed pairs of qubits he knows are described 
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by states pa = {I - Aa)|0)(0| + \a\1){1\,Pb = (1 - 
Ab)|0)(0| + Ai3|l>(l|, < Xa,Xb < \ with Hamilto- 
nians Ha = = and he can makes 

measurements in these energy bases to find the expec- 
tation values (Ha), (Hb)- Charlie knows the states of 
the qubits pA, Pb and therefore also the initial ener- 
gies Ea = TApaHa] = Xa.Eb = Tr[pB77s] = A^. 
What Charlie does not know is the global state p of 
the two systems A and B. Without any knowledge of 
correlations in the system, his description of it is lo- 
cal, that is, Charlie's most unbiased estimate of the 
system state, given his knowledge, is p^*^^ = pa Pb, 
with QMI I{p^'~^'') = 0. He decides to interact A and 
B via an energy conserving unitary interaction, as in 
the earlier heat flow model. If any correlations are 
present, he may see negative heat flow as an indicator 
of them. His description of the system p^ is the mini- 
mally correlated state in its unitary orbit with spectrum 
Aic) = {(1-A^)(1-Ab),(1-Aa)Ab,A^(1-Ab),A^Ab}, 



therefore it must be situated at the minimally correlated 
vertex Rc in its reduced state spectrum set 7^(A('^'), see 
ngure [5] The equal energy region TZe{^''^'') for this state 
is the line that it is situated on, highlighted in the fig- 
ure. All the states corresponding to A^*^)) are obtain- 
able by performing a strong energy conserving unitary 
(SECU) transformation on p^'^', such as f7(°) (9) given m 
equation (19 1. Charlie therefore interacts the qubits so 
that their evolution is governed by U^°\9) and finds that 
he always observes normal heat flow. He concludes that 
his description of the system is correct. 

Now we look from the demon's perspective: it has ac- 
cess to non-local observables and so decribes the system 
via a correlated state p^^^ with the property Tr^(p*^-°)) — 
Tr^(p(^)) = pb and Trs(p(^)) = TrB(p«^)) = Pa- 
To make this example more concrete, we specify the 
state p^^-* with spectrum A^^-* — {Ai, A2, A3, A4}, and 
Ai > A2 > A3 > A4, to be 



/ i(Ai + A4 + (Ai-A4) 



i(A4 — Ai) sin(^ 



cos ( 




A2 








A3 




i(Al + A4 



- Ai) sin 4> 



(Ai — A4) cos ( 



(33) 



and for the demon's reduced states to be compatible 
with Charlie's we demand that ^(1 — A2 + A3 — (Ai — 
A4) cos (f) = \a and 5 (1 -f A2 — A3 — (Ai — A4) cos (j)) — 
Xb- These equalities are satisfied, for instance, when 
A(^' = {0.6,0.3,0.1,0} in which case A^ As - 0.2 



and cos d ~ 



0.3 



last condition requires A^ > 0.1. 



Picking e.g. Xb — 0.4 gives valid A^ = 0.2, cos (/i — |, 
and for this choice of values /(p'^') w 0.4, so the de- 
mon's state is correlated, but has the same marginals 
as Charlie. However under the strong energy conserv- 
ing unitary that Charlie is applying, p^-^^ evolves to the 



state p(Ai, A2, A3, cos 6*, |), see equation (|17|), usmg re- 



sults from section \TV\ one can see that the QMI is con- 
cave over the domain of 9 and its minima occurs at 
the extrema, i.e. when 9 = 0{n), which is p^^K So 
I{U^°^9)p^^'>U'^°\9Y) > I{p^^y)W9 and the heat fiow 
model dictates that normal heat flow will always occur 
for this special correlated state p*^^' under the action of 
U''°\ this is why Charlie always sees normal heat flow, 
even though the system is in reality correlated: his uni- 
tary 1/'^°^ does not "reveal" the correlations. 



However, since p^^'' and p^^^ have different spectra, 
they are located on different unitary orbits and deflne 
different marginal eigenvalue sets 7?,(A^'^^) and 7?,(A^^)). 
Sticking with the values chosen for all the eigenvalues 
above, we see (pictorially in figure [5| that TZ{A^^'>) C 
TZ{A^^^). Furthermore, the equal energy region from the 
demon's point of view is everything one and below the 
thick and thin black line, i.e. 7^i^(A(<^)) C 7^i^(A(^)). 
It is this fact that gives the demon his advantage: be- 
cause of its knowledge of the correlations and the larger 
constant energy region, it can now transform the system 
to a state in 7?,£;(A'^)) that is less correlated and trans- 
fers heat from the cold to the hot qubit. This cannot be 
done using a SECU, instead the demon must use a weak 
energy conserving unitary (WECU) that preserves total 
energy energy on average but does not commute with the 
total Hamiltonian H = Ha + Hb ■ The demon decides to 
interfer with Charlie's experiment and, while Charlie is 
not looking, switches on a demonic unitary interaction, 
producing a state 



/ KAi + A3 + (Ai-A3)cos0') i(A4-A3) 

A2 

A4 

i(A3-Ai)sin(/)' i(Ai + A3-(Ai- 



(34) 



A3) cos ^') J 
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By taking partial traces, one finds that this state has 
thermal marginals (not necessary for our heat flow model 
but included to demonstrate the importance of WE- 
CUs), and spec(7^^^) = A^^^ hence it is on the uni- 
tary orbit of The total energy of 7^^^ is Tt[{Ha + 
Hb)"/''^^] = 1 + (A3 — Ai) cos 0', to satisfy energy conser- 
vation Tr[{HA + HbH'^^] = Ti[{Ha + Hb)p''°^] = 0.6, 
this sets cos(/)' = |. The final QMI is Iij('^'>) « 0.3 < 
/(p*^^)) and system A, which started with energy 0.2, 
now has energy 0.15 and system B, which started with 
energy 0.4, now has energy 0.45. Hence the cold system 
has got colder, the hot system has got hotter and energy 
has remained constant and we have chosen an interaction 
that results in local thermal states so that we can talk 
about temperature (hotness and coldness). 

When Charlie makes his usual energy measurements 
of the state that the demon has tampered with, he ob- 
serves anomalous heat flow and is very confused because 
he thought he had correctly described the system as being 
uncorrelated so only expected normal heat exchange (all 
of his experiments led him to believe this). The observa- 
tion of heat flow from cold to hot with no other change in 
the system (constant energy and constant correlations) 
looks like a real paradox and potentially a violation of 
the second law. This shows that thermodynamics is only 
properly deflned with respect to the particular class of 
observables that are used to prepare and manipulate the 
system. In turn, thermodynamic entropy must be deflned 
with respect to this class if one wants to universally up- 
hold the second law. 



0.5 - 



0.4 - 



0.3 - 



0.2 - 



0.1 ■ 



0.0 - 




0.0 0.1 0.2 03 0.4 0.5 

FIG. 5: The darker (embedded) grey region is the set Charlie 
thinks he has unitary access to 7?.(A^'^^). His description of 
the state is located and the minimally correlated vertex Rc of 
this set, the energy conserving region 7^£;(A'^') for Charlie's 
state is the thick black line. The entire (light plus dark) grey 
region is the set of states the demon has access to, 7?.(A^^^) 
(a different unitary orbit), and its description of the system's 
marginals coincides with Charlie's, but the state the demon 
detects p'^' is now some correlated state in 7?.(A''^'). The 
demon's equal energy region is on and below the thick and 
thin black line. To confuse Charlie it evolves the state p*-^^ 
to 7'^^ which is situated at D. 



VI. DETERMINING p^in IN THE GENERAL 
DIMENSIONAL CASE 

We now return to the problem of determining the abso- 
lute minimum attainable correlation on a given unitary 
orbit for a bipartite system of arbitrary size. For two 
qubits we found Pmin by directly minimising /(p) over 
the unitary orbit. In fact, it is not necessary to solve the 
whole quantum marginal problem in order to find Pmin - 
sufficient is the solution to the simpler problem of finding 
the convex hull of the unitary orbit. 

Even though the domain is now extended to the convex 
hull of Op, we initially ignore the total entropy S{p) in 
our optimisation and focus on just the variation of the 
hmction G{(j) S{(Ta) + S{a'B) over all states a in the 
convex hull. Omitting S{p) is not a problem because as 
it turns out the minimum of G{a) over the convex hull 
will be given by a state Pmin on the unitary orbit of p. 

A state is in the convex hull Cp of Op if it can be written 
in the form a = '^iPiUipU} for a set of unitaries {Ui\ 
and with X^jPi ~ ^^Pi > OVz. An equivalent condition 
is [T5^ spec(cr) =: S ^ spec(p) — A where S,A are 
vectors of the eigenvalues of cr, p respectively and we have 
introduced the majorization relation as a partial order on 
vectors. A vector p is said to majorise another vector q. 



denoted q -< p, for each k = 1, . . . ,d the following holds: 

i2if^i2Pr (35) 

i=l j=l 

d is the dimension of p,q, Pi,qi > and equality holds 
only when k = d [19 . In Ref. [7j it is shown that 
the eigenvalues of all marginal states in Cp are given 
by marginals of probability distributions P majorised 
by A. This reduces the minimisation down to a clas- 
sical problem: the function G{<7) = S{crA) + S{(7b) — 
H{Pa) + H{Pb) -■ G{P), where Pa,Pb are marginals 
of P and are sets containing elements of the vector P, is 
now to be minimised over the convex set of probability 
vectors satisfying P ^ A. It is straight forward to show 
that G{P) is concave on this set of vectors. Therefore 
minima occur at the extrema of this set, that is when 
P = n(A) [7], where 11 is a permutation on the compo- 
nents of A. The states that these extrema correspond to 
are minimally correlated states Pmin G Cp fulfilling the 
property diag(pmi„) — n(A), thus they are in fact on the 
unitary orbit: prmn S Op. 

So we have reduced the problem down to the following: 
the minimum of the QMI on a unitary orbit has a value 
of /™„ = H{[n{A)]A) + H{[U{A)]b) ~ H{A), and the 
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state that corresponds to this is a classically correlated 



B. Minimising the H{Aa) + H{Ab) over 
permutations of the elements of A 



dA.dB 

E ^3.k\j)A{j\(^\k)B{ 
J,k=l 



where {Ij)^}, are local bases and Xj^k is a rein- 

dexing of the eigenvalues A^: Xj,k 



K=U~l)dB+k- 

As in the two qubit case, the ordering of the eigen- 
values matters, we now turn our attention to finding the 
permutation IT that gives the minimum QMI. 



A. The convex hull Cp of the two qubit system 

Before we move on it is constructive to regard what 
the convex hull looks like for the two qubit system. In 
section 



III B we found the minimum of the QMI "manu- 
ally" because the domain TZ was not convex. The convex 
hull of a typical TZ is shown next to TZ in figure [6] for 
< Xa,Xb 1^ \- It is evident now what reduced states 
are absent in the unitary orbit: unitary evolution con- 
strains the reduced states to be closer in mixedness than 
the map that is a convex combination of unitaries. The 
constancy of the joint state spectrum acts as a "tension" 
between the reduced states, restricting the marginals to 
being of similar Bloch vector length. In fact, the inequal- 
ity (Eq. [6]) that cuts away the part of the convex hull so as 
to give TZ is clearly the defining feature for the marginals 
of Op and is the hardest to find. As Bravyi puts it, this 
task is "rather formidable" [7] . It is this lack of convexity 
that contributes to why the quantum marginal problem 
is so difhcult. 



/ 



FIG. 6: Comparing the set of reduced state spectra TZ, on 
the left, and its convex hull, on the right. Although TZ on its 
own looks like a convex polygon, one must remember that in 
fact the eigenvalues Xa, Xb are defined over the range [0, 1], 
whereas here only a portion of this, [0, |], is represented. If 
we imagine (using symmetry) the TZ over the full range [0, 1], 
which is the domain on which the entropies H{Xa), H{Xb) 
are defined, we see that it is not convex. 



From now on, unless explicitly stated, we ignore the 
joint entropy H{A) as we will be looking at individual 
orbits. Let the components of the vector of eigenvalues 
A, which is essentially a bipartite probability distribu- 
tion, be arranged in non-increasing order, > A^+i Vz = 
1, . . . , d — 1. The d — dAds eigenvalues can be arranged 
in a table T, such that the first ds elements of A occupy 
the first row of T, the second ds elements reside in the 
second row, and so on until we occupy dA rows and ds 
columns. This setup makes it easier to see how permu- 
tations of these elements affect the marginal probability 
vectors A^, A^. 

Each eigenvalue A^ is identified with a table entry Tr,c, 
the rows are labelled by r = 1, ... ,0?^ and columns by 
c = 1, . . . , (is, and we denote by R^'') (C'^'^^) the vector 
of table components in the r-th row (c-th column). The 
table below demonstrates this: it shows just one way of 
arranging the A^s, 



r \ c 


1 


2 




ds 


1 


Ai 


A2 




Xds 


2 


Adfj+i 


^dB+2 




^2dB 












dA 
r\c 


1 


A-l) + 

2 


1 


Ad 


B{dA-l)+2 

ds 




^dAdB 


1 




Tl,2 




Tl,dB 


r 


2 


T2,l 


T2,2 




T2.dB = 












dA 


'^dA.l 


TdA.2 




TdA,dB 



and 



C^'^ ^{ri,c,...,rdA,c) 



When written in tabular form, it is easy to calculate 
the marginal probabilities vectors (MPVs) of A. The 
MPVs for systems A and B are denoted a and b respec- 
tively, their components are just the sums of the elements 
in the row and column vectors of the table T 



ds 

1=1 


dB 

1=1 


ds 

■•^E^i' 

1=1 


dA 


dA 


dB 


1=1 


•■>E^1^'' 

1=1 


■•^E^l' 

i=l 



a and b are equivalent to what we previously called 
Aa and A^. We use n(T) to mean the permutation 
on the table elements that corresponds to n(A). A re- 
ordering n changes the MPVs a and b and hence alters 
H (a) + H (b); we seek the permutation 11 that gives the 
minimum of the sum of entropies. 
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1. Sorting rows and columns into decreasing order 

We define two sorting algorithms on the tables: Sf"™ 
(21™') sorts rows (columns) into descending order. Con- 
sider all d\ tables of d elements, which are denoted T^, 
applying these operations in tandem to all tables gen- 
erates a subset Yd of tables with all their elements ar- 
ranged in decreasing order in their rows and columns. 
Yd := {T^ = 2tTd : VTrf,a := 2t™'^2l™'} = {T*^ = 
2l*Td : VT<j,2t* := Sl^'St™""}. Note however that in the 
case when dj^ = ds, we define this set as not contain- 
ing tables that are transposes of another table in Yd even 
though these are generated by the sorting operations, 
that is G ^(^a)^ ^^'^ (T^)'^ ^ Y^dA^- This is because 
they have the same QMI (a transposition on a table is a 
symmetry of the QMI, see section VIC). 



The table that dictates how the eigenvalues should be 
arranged in order to obtain the minimum mutual infor- 
mation is in the set Yd. To see this, consider an unordered 
table Td- The entries of row MPVs a and are 

dB 
i=l 

where R^'^^^ is the r-th row vector of St^'T^. In Ref. [19] 
it is shown that a and a^^ obey a majorization relation a -< 
a.^ and since entropy is Schur concave |19j . this means 
H{a) > (a^). 

In addition sorting the rows and obtaining table StT^; e 
Yd leads also to b ^ b^, and so H (b) > H (b^) leading 
to 

H{a) + H (b) > H (a^) + H (b^) . 

Hence the operation 2t (or 2t*) lowers or keeps constant 
the sum of entropies and the minimum of the mutual 
information is always given by a table in Yd- Note that 
this also holds if 21''°"' , 21"^°' are defined as putting row 
and column elements into ascending order. 

The set of tables that result from this sorting are the 
dA ^ dB Young tableaux. The number of these tables is 
given by the hook formula of Frame- Thrall- Robinson [20] 



C. An Aside: Symmetry of the Mutual 
Information in Terms of Tables 

We noted earlier in section [Tij that the quantum mu- 
tual information is invariant under particular transforma- 
tions, they were local unitary transformations and a swap 
of subsystem states. We can now talk about this in the 
language of table element permutations. The notation 
/(T) symbolizes the mutual information of the probabil- 
ity distribution arranged in the order described by table 
T. Transformations on the table T which are symmetries 
of the QMI are: 

1) Transposition (only for tables where dA — ds) 

I (T^) = / (T) 

This amounts to swapping the distributions of A and B. 

2) Permuting rows: writing the table with row vectors 
as T [R*^^-', . . . , R^'''^''] , and a permutation on the row 
vector ordering IIr 



/(T) =1 (t R(i\...,R('*'*) 
Hk (r(i ),..., Rt''- 



/ T 



This is a local operation on the subsystem B. 

3) Permuting columns: writing the table with column 
vectors as T [C(^\ . . . , C'^'*^^] , and a permutation on the 
column vector ordering He 



c(i),...,c('^«^ 



/(T) = / T 



f~,{dE 



This is a local operation on the subsystem B. 

4) Swapping all table components r^s from, say, de- 
scending to ascending order: a combination of the above 
transformations. 

These transformations can be classed in terms of a 
group action. The symmetry group for / (p) is the per- 
mutation group Sd oi order D\ acting independently on 
the rows and columns. In a d-dimensional system there 
are d\ permutations/tables. The symmetries reduce the 
number of tables by a factor of dA^-dsU and an additional 
factor of 2 if = d_B . Note that the Young tableaux are 
a further subset of this reduced set. This is discussed in 
more detail in Ref. [3]- 



N{d, 



n 



j=dA 



(36) 



D. Special Case: Minimum mutual information 
state when dA = rfs = 2 



and this series grows exponentially: N{2, 2) = 2, 
N{2, 3) = 5, iV(3, 3) = 42, iV(4, 4) = 24024. 

Our set of interest is Yd, and because it is defined as 
not containing transposes in the case when dA — ds, 
there are in fact only ^N{dA,dA) tables in Y(^dA)^' these 
are what we call the "independent Young tableaux" . 



For a two qubit system, there is one independent 
Young table, from Eq. ([36|, ^^^(2, 2) = 1, it is 



r \ c 


1 


2 


1 


Ai 


A2 


2 


A3 


A4 
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This table minimises -ff (a) + H{h) for all A and corre- 
sponds to the state found in equation ([9]), confirming the 
derivation in section UlIBI 

In total there are 4! — 24 tables of this dimensionality, 
but due to symmetry of the QMI there are only 24/2^ — 3 
independent tables, i.e. only three distinct values of the 
QMI. The tables chosen to represent them are the one 
given above and 



T 



(1) 



r \ c 


1 


2 


1 


A2 


Ai 


2 


A3 


A4 



r \ c 


1 


2 


1 


A3 


Ai 


2 


A2 


A4 



Note these two tables are not in Y4. 

We can order the row and column MPVs of these tables 
with respect to majorization 



(1)^ - 



(37) 

L4 ) — LFv^4 ) ~^ "^^4) (38) 

Writing the QMI of table T as /(T) we have then that 



b(Tf ') -< b(Ti) 



(39) 



This gives us a full ranking of the mutual information in 
the full two qubit quantum space. On a unitary orbit, in 
increasing order 

/(mm, class) = H{\3 + A4) + i? (A2 + A4) - H{A) 



I{max, class) = H{X2 + A4) + H{Xi + A4) - H{A) 



some A. For instance, when cIa = 2, = 3, the dis- 
tribution A = ^{6,5,4,3,2,1} minimises the MI when 



(3) 



the elements are arranged as in table Tg 
A = ^{10,9,8,3,2,1}, then the optimal table is T, 



whereas if 
(1) 

6 ' 



these tables are given in section VI E 1 

Below we list explicitly the Young tables for the cases 
when = 2, = 3 and cLa ^ Ab = 3. The set Yd 
is the collection of "independent" (not including trans- 
poses) Young tables T^, « = 1, . . . \Yd\ — N . From now 
we use the less cumbersome notation i := A,, i = 1, . . . , d. 
Also included are frequency plots in figure[7]for the tables 
that minimise the MI for randomly generated probability 
distributions. 




1 1.5 2 2.5 3 3.S 4 4.6 S / 



limax, Sep) < logd^i — 1 



I{max, quant) — 2 — H{A) 

where the QMIs are over the set of states that are class = 
classical, sep — separable, quant = arbitrary quantum. 
The maximum separable value is ascertained from the 
upper bound of the classical mutual information over the 
set of all probability distributions with four elements and 
assumes 2 — H{A) > 1, or equivalently H{A) < 1. This 
full ordering of the mutual information is possible only 
for two qubits. 



E. The correlations structure in higher dimensions: 

dA >2,dB>3 



X 10 




20 I 



For systems of larger dimension, the situation becomes 
more complex because there are now multiple Young 
tableaux from which we must determine the one that 
is the minimiser for the mutual information (MI) (the 
problem is now purely classical). 

A general solution for d > 4 does not yet exist and, 
unlike the 2x2 case, it depends on the choice of A itself; 
every member of the Young tableaux set is optimal for 



FIG. 7: Histograms displaying the number of times c a ta- 
ble T^'l^^^^ (labelled i) minimises the mutual information 
for d — 6 top figure and d = 9 bottom figure. In each case 
1,000,000 full rank probability distributions were generated 
randomly. Surprisingly the most frequently occurring tables 



are Tg'*' and Tg 
togram. 



(19) 



in the latter case Tg^^' dominates the his 
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1. The distinguished set of tables for dA = 2, d_B = 3 : Ye 



For the qubit-qutrit direct calculation results 

in the following set of tables for Yg 



1P(1) 

^6 



1 


2 


3 


4 


5 


6 



■ Tr(2) _ 



1 


2 


4 


3 


5 


6 



: ^(3) _ 



1 


2 


5 


3 


4 


6 



Jig — ■ 



1 


3 


5 


2 


4 


6 



Tr(5) _ 

Jig — ■ 



1 


3 


4 


2 


5 


6 



Curiously, although we cannot deduce which table in 
Yq gives the minimum MI for a given A, we do know 
there is a unique table (not in the set Fg) that gives the 
maximum classical MI. The table is 



Tnax, class 



1 


4 


5 


6 


3 


2 



T 



and this is proved in Ref. [3 by introducing a new kind 
of partial order called "quasi-majorisation" . 

So for cIa — 2,dB — 3, we only partially know the 
ranking of the mutual information: 

I{min, class) < I{max, class) 



I{max, class) = H{X2+\4+Xe)+H{Xi+X6, X3+\4)-H{A) 



T, 



(1) 



T, 



(4) _ 



T, 
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T 
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T 



(19) 
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For this, and larger dimensions, the minimum (and 
maximum) classical mutual information values are not 
known, there is not a unique table that defines these be- 
cause it depends on the spectrum A. 



A Partial Ordering on the Young Tableaux and 
the "See-Saw" Effect 



I{max, class) < I{max, sep) < log(min{(i^, d^}) = 1 



I{max, quant) < logG. 

Broadly speaking, it is possible for a minimally corre- 
lated classical state to be unitarily evolved to the max- 
imally correlated classical state, however this classical 
state is not necessarily the most correlated state over 
the set of separable states intersecting the unitary orbit. 
Thus a further unitary can increase the correlations to 
form a separable state which has non-zero discord. Be- 
yond this, it may be possible to leave the convex set of 
separable states and obtain an entangled mixed state. 



2. The distinguished set of tables for dA = ds = S : Yq 

For the case of two qutrits we find that the set Yg is 
given by 



We cannot say whether the sum H (a) + H (b) will 
increase from one table to another, nevertheless it is pos- 
sible in some cases to determine how the individual en- 
tropies will change. In the following, we assume a priori 
majorizatioii relations between vectors, that is partial or- 
dering in the general case without reference to a specific 
probability distribution A. 

(X) 

Lemma If two cLa ^ = d Young tableaux , 

( 7) 

T)^ e Yd are related by one transposition 

Tf -^(rs:t^.)Tf 

where tt [rs : tu) swaps the elements t^s and Tm of table 
T^;'^'' , then either 

if (5 > or 

a(^)^a(-^),b(^)^b(-^) 

if (5 < , where a^'^-' and b^'^-' are the row and column 
MPVs of table T^-^) ,X€{T,J), and <5 = Tt„ - r^^. Say 
r > t in the following, without loss of generality. 
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Proof The swap operator tt (rs : tu) cannot arbitrarily 
swap any two nonequal elements of a Young tableau and 
stay within the set Yj^. For a valid swap, it is necessary 
that 



r > t 
s < u 



(40) 



If either r = t or s = u, then the swapping occurs between 
elements in the same row or column, which is prohibited 
because it would violate the Young tableau element or- 
dering, similarly ii r > t and s > t. Note, however, that 
not all swaps which satisfy equation (40 ) are valid swaps, 
for example when cIa — ds — 3, swapping ri2 and T21 in 
table concurs with these conditions but generates a 
table that is not in Young tableaux form. Therefore con- 
ditions (40) are necessary but not sufficient. 



Consider tables T~, 



TT (rs : tu) ¥P 



Two components of the row and column MPVs of 
are affected by the swap, they are 



Yd. 





= Trl + ■ 


■ ■ + Trs + ■ 


■ ■ + TrdB 


(I) 
at 


= Tti + ■ 


■ + Ttu + ■ 


■ + TtdB 


6f) 




■ + Trs + ■ 


■ + ^dAS 


5(x) 


= Tul + - 


■ + Ttu + ■ 


■ + TdAU 



After the swap, the MPVs of 
terms of a*^-^' and b^-^' components 



(J) 



can be written in 





= Trl + ■ 


■ + Ttu + ■ 


■ ■ + TrdB 


= af) 


+ S 




= Ttl+- 


■ + Trs + ■ 


■ + 'TtdB 


. af) - 


-S 




= Tls + ■ 


■ +Ttu + ■ 


■+TdAS 




+ 6 




= Tlu + ■ 


■ ■ + Trs + ■ 


■ ■ + TdAU 




-5 



('T\ (T\ 

where S = Ttu — Trs and \S\ < a\ ,bj and i = 1, . . . ,dA, 
j = 1, . . . , de- The latter condition on the magnitude of 
6 ensures that adding and subtracting it to the compo- 
nents of the MPVs preserves their order; this is important 
when comparing the majorization relation between two 
vectors. If the elements of a^-''^ and b^^-^^ are ordered in 
decreasing order, then the elements of a^-^' and b''-'^' will 
be also. This is a consequence of the fact that we only 
consider swaps which are operations f Yd Yd and 
more importantly that the arrangement of elements in 
the Young tables ensure the MPV elements are ordered; 
for row MPVs ai > 02 > . . . ad a and similarly for column 
MPVs. 

Collating the aformentioned facts, it can be shown that 
the MPVs satisfy the majorization relations 



if (5 > and 



if (5 < 0. □ 



a(-^) ^a(^),b(-^) ^b(^) 



Since entropy is Schur concave, this means that for, 
say, 6 > 



Ha 



> H 



(a(^)),i/(b(-^)) <i/(b(^)) 



and because a transposition causes one entropy to in- 
crease and the other to decrease in this way, we call this 
the "see-saw" effect. 

It is possible to represent this effect by a directed graph 
G = G{Yd). The nodes correspond to the tables and 
the edges connect two tables that are related by a swap. 
A graph G can only represent row or column MPV ma- 
jorization; we call the two types of graphs Grow and Gcoi- 
The direction of the edge connecting two nodes represents 
an MPV majorization relation, hence an arrow "— from 
node I to on the graph Grow implies a'-^^ >~ a^-^K 
From the lemma above this also means that b(^' ^ b(-^) 
is also true. Hence on Grow, a "— arrow will exist be- 
tween nodes I and J^, whereas on Gcoi there will be a 
"■s— " arrow between the same two nodes. Hence Grow 
and Gcoi are equal up to a change in their arrow direc- 
tions, this allows us to just focus on either row or column 
majorization since one graph is easily obtained from the 
other. Figures [S] and [9] display the graphs Grow {Yq) and 
Grow (Yg) ■ 




Node 3 




Node 5 


\ 






FIG. 8: Connected and directed graph Grow (Ya), ea ch nod e 



represent a table, e.g. Node 1 is table Tg^' in section 



VIEl 



Because of the transitivity property of majorization, if 
for tables T^^) , T^-^), T^^) € Yd we have that a(^) >■ a^-^) 
and a(-^) y a^^), then a(^) >- a^^). The graph Grow {Yd) 
captures this more general majorization relation. Follow- 
ing a sequence of nodes connected by edges all directed 
the same way gives a subset of tables that obey this par- 
tial order. For example, in the graph Grow (Yq) , a set 
of connected nodes is Node 1 Node 2 — >• Node 3 — > 
Node 4 from which we can say, for instance, that a^ >- a^ 
and ^ b*. Except for Nodes 3 and 5, all other pairs 
of nodes in Grow [Yq) are connected by arrows all point- 
ing in the same direction hence we can determine almost 
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|Nocle12| 


1 Node 20 1 








Inoa 


B14| 




= 17| 



[Mode 1 9 I 

FIG. 9: Connected and directed graph Grom (Yg), ea ch node 



represent a table, e.g. Node 1 is table Tg'^' in section 
It would be symmetrical like Grow (Ya) if the transposes of 
the tables in Yg were included. 



VIE2 



all the column and row MPV majorization relations of 
any two tables in Yg. Upon inspection of the tables cor- 
responding to Nodes 3 and 5, we see their MPVs are 
incomparable, which is reflected by the fact they are un- 
connected in G(Yq). Hence we can make the following 
statement: 

Two tables T*^-^^ and T'--'^^ in Yd, represented by nodes 
N^^^i and N^'^^ in the graph Grow {Yd) , riave a majoriza- 
tion relation between their row and column MPVs if N'-^^ 
and N^^'' are connected by arrows all pointing in the 
same direction. 



G. Special Case Results 

In some cases we can use simple mathematics to cal- 
culate whether the mutual information will increase or 
decrease from one permutation to another. Let us denote 
a probability vector with n elements p = {pi,p2, ■ ■ ■ ,Pn) 
where pi > P2 > • • • > • Two results follow. 

Lemma 1 Given 
T(i) 



Pi 




Pk 


Pk+i 










Pn 



Pi 




pk+1 


Pk 










Pn 



with fc < n, and p^ = for all i > k. Then /(tI,^^) < 

/(tI"'). 



Proof 1 Let h{x) = —xlogx. Then (ignoring the total 
entropy contribution to the mutual information) 



i{r'i^)^Y.h{p^) 



(41) 



fe-i 



^k-l 



i=2 \i=l / 

(42) 

hence it must be shown that 

/(tP)-/(tW)>0 (43) 



or 



h{pi+Pk) + h{l-pk)-h{pi)>{) (44) 

Writing this as h{pi +pk) — h{pi) +/i(l —pk) — h{l) and 
using the mean value theorem fib) — f{a) = (b — a)f'{c) 
where c e [a, b] then the condition becomes 



/i'(pi-f ^)-/i'(l-a) >0 



(45) 



where a,/3 G [Q,pk\- The derivative of the function is 
h' {x) — — fcilna; — /c2 where fci,fc2 are constants so we 
obtain 



In > 

Pi+P~ 



I~a>pi+I3 



(46) 



(47) 



In the worst case (minimum RHS and maximum LHS) 
the inequality becomes 



1 - Pfc > Pi 

which is true since 1 = ^ Pi +Pfc- D 

Lemma 2 If we have 

T(i) = 
t(2) ^ 



(48) 



Pi 




Pk 


pfc+1 










Pn 



Pi 




Pk + 1 


Pk 










Pn 



with k <n, and p, = for alH > fc + 1. Then L{TiP) < 
/(Ti'^). 



Proof 2 The condition 

/(T(f))-/(T«)>0 



(49) 
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is equivalent to 

Hpi +Pk) - Hpi +Pk+i) + h{l -pk) - h{l -pk+i) > 

(50) 

Using the mean value theorem again we obtain 

h'{pi + Pk+i +l)-h'{l- pk + 6)>0 (51) 
where e [0,pk -Pfe+i], so 

In + ^ >0 (52) 

Pi + Pk+i + 7 

or 

1 -J3fe + <5 >Pi +P/C+1 +7 (53) 

In the worst case (minimum RHS and maximum LHS) 
the inequality becomes 

1 - Pfc >Pi+Pk (54) 

which is true since 1 — J2i=i Pi ^ Pi + Pk-i + Pk > 
Pi + 2pk. □ 

VII. CONCLUSION 

Intuitively it would seem right that the state on a uni- 
tary orbit that is the maximally (minimally) correlated is 
the one which is a convex combination of maximum (min- 
imum) QMI pure states. Here we have formally shown 
that this is true, however it is highly non-trivial and addi- 
tional considerations beyond majorization are required. 
The entire analysis boils down to the dependence of the 
mutual information on orderings of elements in Young 
tableaux, and these tables display a subtle property in 
terms of a partial order with respect to marginal en- 
tropies. 

This work finds application to closed system thermo- 
dynamics when one includes the additional constraint of 
energy conservation. 

Many challenging questions remain, the obvious one 
being, given a spectrum of a bipartite system, what is 
the arrangement of its eigenvalues that gives the mini- 
mum (and maximum classical) mutual information? Can 
anything be gleaned from our analysis about other cor- 
relation/entanglement/concave functions? What about 
the inverse of the problem, varying the QMI when the 
marginals are fixed but the global state is not (but must 
be compatible with the local description)? 

In terms of thermodynamical work, can we devise a 
simple collision model in which correlations can remain 
non-zero but equilibration still occurs? Can our work be 
connected to recent advancements in the resource theory 
of athermal states I21j? What do our results contribute 
to traditional open systems quantum thermodynamics? 

We see that there are many avenues that our investiga- 
tions could lead on to, and hope that our work provides 
use to people investigating the thermodynamics of quan- 
tum systems. 
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Appendix A: Details on the construction of marginal 
spectras 

For two qubits, the minimum QMI state, Pmin, must 
correspond to a point on the line (jp, p = 1, 2 in figure [T] 
There are two cases to consider. 



1. p = 1, A2 = A3 

The inequality in equation ([s]) reduces to 

+ As > 2(A3 + A4) 

and it is now redundant because it is no longer indepen- 
dent - it can be formed by adding inequalities ([3| and 
Q. Therefore it places no additional constraint on TZ 
and the boundary closest to the origin is in fact a cor- 
ner where the lines /i and /2 intersect. The minimum 
values of \a and As occur at this intersection, where 
Aa = As = A3 + A4. 



a. p — 2, A2 > A3 

Now 

Xa + Xb> 2(A3 + A4) 

and this inequality is a true constraint on TZ. The line 
92 cuts the bottom left corner of TZ and the minimum of 
H{Xa) + H{Xb) is somewhere along it. Let us define 

K = \2 + X3 + 2A4 

then g2 can be rewritten Xb — K — Xa and the function 
to be minimised on this line is H{Xa) + H{K — Xa), with 
Xa € [0, K]. The domain is now convex and the function 
concave and symmetric about the maximum at = ^ 
(where Xa = Xb, the symmetry line of TZ) hence the 
minima occur at the edges of 172 ■ 

Let us call the coordinates of the boundaries of 52 
-Ri,i?2 G T^, shown in figure [ij By symmetry we need 
only consider one of these points, Ri say (i?2 is given by 
swapping Xa,Xb in -Ri), its coordinate is either at the 
intersection of the lines hi and (72 or /2 and 32- /2 and 
52 meet at {Xa,Xb) = (A2 -I- A4, A3 -I- A4) =: (A^, A^). 
However there are two situations to consider for the inter- 
section of g2 and hi depending on min{Ai — A3, A2 — A4}. 
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If hi = miii{Ai — A3, A2 — A4} = Ai — A3: 
then 52 and hi meet at (A^, As) 

= (^(Ai - A2) + A2 + A4, A3 + A4 - ^(Ai - A2)) 
=: (A^, X'q). 



As noted above, the point Ri in TZ corresponding to the 
minimum QMI state is an intersection of two boundary 
hues: A^i = A3 + A4 and A^ + A^ = A2 + A3 + 2A4. Hence 
Pmin must commutc with the two operators O, O' that 
give rise to these hnes. 

To obtain inequahty Q we set O = Ha, then 
{01,02,03,04} = {1,1, 0,0} and 



If hi = min{Ai — A3, A2 — A4} = A2 — A4: 
then g2 and hi meet at (Aa,A_b) 



- (2(A3-A4) + A2 + A4,A3 + A4 

(A:;,A'^). 



;(A3-A4)) (Al) 
(A2) 



It is clear that A'4,A4 > A*i and 



< At 



This imphes that in both cases hi intersects g2 at a 
point that is further from the diagonal Xa = Xb than 
the intersection of /2 and 52 and so the coordinates 
(A;4, A^), (A'^, A^) are not in 7^, hence Ri = {X*^,X*b). 

Therefore, the minimum of H{Xa) + H{Xb) occurs at 
Ri where 



Xa = tHp^A) 

> inf tr(rn^) 

= A4 + A3 



(AlO) 
(All) 

(A12) 



Similarly, to obtain inequality ([5| we set O' — Ha + 
Hb- Let us suppose, without loss of generality, that Ha — 
|0)(0| and Ub = \v){v\, where \v) = Ub\0) and Ub is a 
single qubit unitary. O and UgOUB are cospectral and 
have eigenvalues {Oi, O2, O3, O4} = {2, 1, 1, 0}. Thus 



AA + Ai 



= tr(p(nA 4 
> inf 

re-D('H,4.A) 



tr{T(nA 



2A4 + A3 + A2 



(A13) 
(A14) 

(A15) 



and at i?2 



Xa — X2 + A4 
As = A3 + A4 



(A3) 

Without loss of generality O, O' can be set to be di- 
(A4) agonal in the same, say computational, basis. In matrix 
form they are 



Aa — A3 + A4 
As = A2 + A4 



(A5) 
(A6) 



It is evident now that gi in the previous is a special case 
of this when A2 = A3. 

The only task left is to determine the composite state 
Pmin which has reduced states with the eigenvalues given 
in equations ( |A3[) and (A4) (or (A5) and (A6)). The in- 
equalities (|3| - ([5|) which give the lines fp,gp are obtained 
by using the relation 



inf tr(OT) 

re-D(-H,4,A) 



A1O4 + A2O3 + A3O2 + A4O1 (A7) 



where O is an arbitrary two-qubit hermitian operator 
with eigenvalues Oi > O2 > O3 > O4 and r is four 
dimensional state with eigenvalues A = {Xi}j^i ordered 
in the same way. This infimum is reached when O and 
T are diagonal in the same basis with their eigenvalues 
arranged in such a way that equation (A7) is obtained. 



Let us write the marginals of a two qubit state p as 

PA = XAnA + {i-XA){i-nA) (A8) 

PB = XbIIb + {1- XB)iI -TIb) (A9) 



o = 



O' 



1 0' 
10 

0, 

2 
10 
10 




(A16) 



(A17) 



Then pmin must satisfy [p. 
0. Alternatively, this 



n,0'] = 

require 

, O + O'] = and O 7^ -O' (which is true). Given 



min,0] = and [p, 
fulfilled if 



IS 



we 



+ 0' 



/ 3 

2 

10 

Vo 



(A18) 



then Pmin must be diagonal in the same, computational, 
basis. Since tr(p,„i„nA) = A3 + A4 then 



/Al 

A2 

A3 

V A4 



(A19) 



and 11^,113 are the projectors onto the eigenstates of 
PAtPb associated with the lowest eigenvalue. 



This state and its equivalence class is a classically cor- 
related state, and it is the minimiser for the QMI. 
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Appendix B: Maximum A/b 



To see more clearly that equation ( 15 ) describes a line 
let us define Xa — y, Xb — x, cosOa = t, cos 0b = s and 
C = 1 — E > 0. Then the equation can be recast in terms 
of these variables 

yt + XS ^ ^{S + t) - C 

We call this line F{s, t). A special case occurs for F(0, 0) 
where E = 1 and no other energy is allowed. These lines 
have symmetry about the lines x — — \- 

The region of states with initial and final energies equal 
to E will be where the lines F{s,t) intersect the region 
TZ. There are four important lines that frame the set of 
all F{s,t) and they occur when s = ztl,t = ±1. The 
intersection of these lines, denoted by F{s,t) F{s' ,t'), 
is given below in terms of the coordinate representation 
{x,y) 



This is depicted in figure 10 A general line in F{s,t) 
passes through the points (Y, 5 ~ t) ^^'^ (5 ~ ?' 5) 
and so it never goes inside the region bounded by the 
lines F{1, 1), F{1,-1), F{-1, 1),F(-1,-1). This is also 
shown in the figure. 



i^(l,l)&^^(l,-l) 



i^(l,l)&i^(-l,l) 



F(-l,-l)&i^(l,-l) 



F(-l,-l)&i^(-l,l) 



1 1 

2' 2 
1 
2 



(Bl) 
(B2) 
(B3) 
(B4) 



The region TZ is restricted to x,y < ^. The set of lines 
F{s, t) that fall into this range are F{1, 1) and all the lines 
"below" it, that is all the lines satisfying x + y < 1 — C. 
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